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A bstract 


The light-like cusp anomalous dimension is a, universal function in the anal- 
ysis of infrared divergences. In maximally UV — 4) supersymmetric Yang- 
Mills theory (SYM) in the planar limit, it is knovn, in principle, to all loop 
orders. The non-planar corrections are not knovn in any theory, vrith the 
first appearing at the four-loop order. The simplest quantity vvhich contains 
this correction is the four-loop tvvo-point form factor of the stress tensor mul- 
tiplet. This form factor vvas largely obtained in integrand form in a previous 
vvork for /V — 4 SYM, up to a free parameter. In this vvork, a reduction 
of the appearing integrals obtained by solving integration-by-parts (TBP) 
identities using a. modifted version of Reduze is reported. "The form factor 
is shovrn to be independent of the remaining parameter at integrand level 
due to an intricate pattern of cancellations after IBP reduction. Moreover, 
tyvo Of the integral topologies vanish after reduction. The appearing master 
integrals are cross-checked using independent algebraic-geometry techniques 
explored in the Mint package. The latter results provide the basis of mas- 
ter integrals applicable to generic form factors, including those in Quantum 
Chromodynamics. Discrepancies betvveen explicitly solving the TIBP rela- 
tions and the MINT approach are highlighted. Remaining bottlenecks to 
completing the computation of the four-loop non-planar cusp anomalous 
dimension in /V/ — 4 SYM and beyond are identified. 
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1. Introduction 


Maximally (UV — 4) supersymmetric Yang-Mills theory (SYM) offers 
perhaps the best chance of truly solving an interacting four-dimensional 
quantum field theory. In addition, it is a proven stepping stone to pio- 
neer nevv computational technology, as the large amount of supersymmetry 
renders many computations much simpler than their non-supersymmetric 
counterparts. In many cases, the resulting technology has transformed com- 
putational povver in generic quantum field theories. For instance, in QCD, 
recent years have seen a, boom of next-to-leading-order (NLO) computa- 
tions considered unfeasible only a, fevv years before, inspired ultimately by 
VVitten s tvvistor string İli). 

The focus of this paper is the computation of the Sudakov form factor, 
vrhich is an observable that involves one gauge-invariant operator from the 
stress tensor multiplet and tvro on-shell massless states, 


7 — (pi, pəlOl0) . (1) 


m.V — 4 SYM, this form factor vvas first discussed and computed at the 
tvvo-loop order in İ2). It is notevvorthy that the three-loop form factor vvas 
first computed in QCD İ3), almost three years after the master integrals 
vvere found in /4İ. "Fhis result vvas then adapted to provide the three-loop 
ansvver in /V — 4 SYM in bl. The latter computation cleanly shovvs that 
the computation vvithin /V — 4 SYM is technically much more straightfor- 
vvard, especially vrhen employing modern unitarity-based methods instead 
of Feynman graph techniques. 

Since there is only a single scale in the problem of computing the Sudakov 
form factor, the dependence on this scale is fixed by dimensional analysis. 
Hence, the form factor evaluates to a, function vvhich only depends on the 
dimensional regularization parameter € (D £ 4 — 2e), the coupling constant 
and the number of colors /V.. Moreover, the dependence on € is governed 
to a large extent by the knovrn structure of infrared (IR) divergences. The 
divergent structure through R.a for instance, is determined by the so-called 
cusp (or soft) anomalous dimension, vvhich gets its name from its appearance 
in the computation of the light-like cusped VVilson line İ6, Hr. The cusp 
anomalous dimension is a, function vvhich is universal for a, given quantum 
field theory. 


m /V — 4 SYM, an integral equation vvhich determines the leading pla- 
nar part of the cusp anomalous dimension in principle exactly vvas derived 
from integrability “Sİ. This is to date the most povverful and precise mani- 
festation of the cluster of ideas knovrmn as the AdS/CFT correspondence İ9). 
Hovvever, much less is knovm about the non-planar part of the cusp anoma- 
lous dimension. At vveak coupling, an immediate problem is that the first 
non-planar correction appears at four loops and has never been calculated 
to date, in any theory. PFurther motivation comes from a, con)ecture 10 
that this correction may in fact vanish. This is based on the assumed com- 
pleteness of the dipole contributions, vvhich is consistent vvith the analysis 
of IR divergences in Yang-Mills theories explored through three-loop order, 
see also . Hovvever, there is evidence from the Regge limit analysis in 

that this naive dipole summation is not sufficient at four loops. İt is 
important to settle this issue by explicit computation. 

The integrand of the form factor vvas obtained in using color-kinematic 
duality İ4, hi he) as an ansatz generator. "he coefflicients in the ansatz 
vvere fixed by a choice of unitarity cuts, up to one remaining parameter. "To 
evaluate the integral, a mayor next step is to perform the integral reduction, 
in particular using integration-by-parts (TBP) reduction İz, . This turns 
out to be highly non-trivial due to the complexity of four-loop integrals. "The 
success of this reduction is a, mafor achievement of this paper. İn particular, 
this reduces the challenging non-planar cusp anomalous dimension problem 
dovvn to the computation of an explicit set of master integrals. 

Based on the TBP reduction result, surprisingly, the free parameter in 
the integrand turns out to be truly free: it drops out of the full form factor 
after the reduction. "his completes the determination of the integrand of 
the /V — 4 SYM form factor at four loops. İn other vvords, vve prove that 
there is a, one-parameter family of four-loop form factor integrands vvhich 
satisfy color-kinematic duality. 

The master integrals are also explored independently by algebraic tech- 
niques. VVhile the results are consistent vvith the IBP reduction in malority, 
some rare examples of mismatch present interesting open questions to under- 
stand further. The results from algebraic methods, vrith some qualifications, 
provide the master integrals of four-loop form factors in generic quantum 
field theory, including QCD. 

The remaining problem of computing the four-loop form factor in /V — 4 
SYM is the evaluation of the master integrals. Many techniques are available 
to compute Feynman integrals, see and references cited therein for an 
overvievv. Hovvever, the basis integrals of the four-loop form factor are rather 
non-trivial. They include non-planar four-loop integrals vvith up to quadratic 


irreducible numerators, the ma)ority of vvhich, to our knovvledge, have never 
been integrated. VVe present partial numerical results and briefly survey the 
challenges involved in completing the computation, vrhile leaving the full 
evaluation to future vvork. 

"This article is structured as follovs. First, the connection betvveen form 
factor and cusp anomalous dimension is revievved in section This sec- 
tion also introduces some of the language to study Feynman integrals and 
discusses the issue of graph symmetry. 1n section İ3) vve report on the IBP 
reduction of the Feynman integrals in the four-loop form factor in /V — 4 
SYM based on the Reduze package İ0). The independence of the sum on 
the free parameter reported in is elucidated. VVe then proceed, in section 
İl to study the basis of master integrals using algebraic techniques explored 
in . A brief survey of the numerical evaluation is given in section İ5l Ap- 
pendices contain information on a three-loop cross-check and an alternative 
choice of basis integrals. ln the supplementary material section, there are 
several ancillary files vvhich contain results for the basis of master integrals. 


2. Reviev and setup 


In the folloving, extensive use vvill be made of the results in for the 
four-loop form factor in /V — 4 SYM, in particular of the graphics and tables 
in section 5 of that paper. The numbering of the integral topologies refers 
to that paper. 

ln the course of vrriting this paper, a typo vvas discovered in Table 5 
of Hs): the entry for the color factor of integral topology 26 should be 
24N2öa,aş. In particular, this integral topology does not contribute to the 
planar cusp anomalous dimension. 


2.1. Relation betuveen form, factor and cusp anomalous dimension 
The Sudakov form factor is an observable vvhich here involves one gauge- 
invariant operator from the stress tensor multiplet and tvvo on-shell massless 
states, 
72 — (pı, pəlOl0). (2) 
In the present paper, vve focus on the four-loop Sudakov form factor in /V — 4 
SYM, for vhich the integrand vvas obtained in 


Hm A — 4 SYM, form factors vvere first studied thirty years ago in and revived 
in the past fev years at vveak coupling İB, bə, bs, ba, bzi, İd bi psi, lg, Bil and 


at strong coupling ERİ 34), 35). See also recent developments concerning form factors of 
non-Bogomolnyi-Prasad-Sommerfield (BPS) operators İd, İB7, ERİERB 


The Sudakov form factor for the stress tensor multiplet can be vritten 
as the tree-level form factor times a scalar function, 


7, — FÜ F.(g, N.). (3) 


It is this function F2 vvhich vill be computed. In momentum space, F3 is 
a function of the only scale in the problem, q? — (pı -- pə)”, the coupling 
constants, group theory Casimirs and e. Here, q is the momentum associated 
vrith the gauge-invariant operator. Note that it does not obey a, mass-shell 
condition, in contrast to the momenta of the gluons, p? — pş — Ü. 

The cusp anomalous dimension “eusp iS related to the tvvo-point . 
factor through the universal exponentiation of IR divergences mim mimimi 
vrhich, in dimensional regularization, takes the form 


0) 
LogİFzl — yeyə” ə m (€). (4) 


The cusp anomalous dimension can be calculated as a, perturbative series in 
the coupling constant. Displaying the dependence on group theory quanti- 
ties up to fourth order, it reads 


ÇI "/cusp ”” . gönün m aıg”CA—aşgtCa Laşg”CA-g" (ağ 62 -. aq”da)-O(g “7” 


l 
(5) 
vrhere CA is the Casimir vvhich, for the SU(.N.) gauge group, simply has value 
N,. The quantity daq is a particular group theory Casimir invariant vvhich, 
R for SU(Nz), .—. the value N£ F 36N2. In previous vrorks, the imimim: 
.. az, aş and al vvere calculated in perturbation inv lid 
r ə and using integrability İsl. m ”t Hoofts planar limit 51 N, ". Co r 
g?N: fixed, it is clear that the first subleading color correction to the cusp 
anomalous dimension occurs at four loops. 

From the general exponentiation formula in equation (4), this non-planar 
correction can be isolated from the non-planar part of the form-factor at 
four-loop order at order R. in the € expansion. Since the integrals present 
in the form factor naively diverge as - this gives an important check on 


2"The three-loop cusp anomalous dimension in /V — 4 SYM vvas first obtained (con)ec- 
turally) in by extracting the leading transcendentality part of the QCD result bi bd 
m this case, the cusp anomalous dimension is related to the anomalous dimension of the 
tvvist-tvvo operator in the limit of spin / — oo. 


any computation of the non-planar form factor: the first six orders in the 1 
expansion must vanish. 

The planar correction to the cusp anomalous dimension at four loops 
can be expressed as a sum of rational coefficients times either ((3)7 or z6 
las, lad, ir, vrhich are both constants of franscendentality six. Based on the 
maximal-transcendentality principle 48 , it is natural to conğecture that the 
four-loop non-planar correction aN)” is given by a combination of the same 
transcendental numbers, if not vanishing. 


2.2. Feynman integrals 


An 12-loop Feynman integral vvith nz so-called indices aı,..., an 1S an 
integral vvhich can be vrritten as 


HR Fe” ..dDiz ((/Diy" ...(U/DA)", (6) 


vrhere İD), are inverse propagators. For future reference, vve define the ab- 
solute value of the sum over positive and negative indices to be £ and s 
respectively. The parameter s is the numerator povver and £ is the propaga- 
tor povver. 1n this article, the indices are strictly integer valued. Integrals 
vvhich have the same non-zero indices for propagators are said to belong to 
the same sector. ln the special case vrhere all propagator povvers are one, 
the integral is referred to as the corner iöntegral of the sector. 

The problem under study in this article consists of a, series of four-loop 
tvvelve-propagator integrals. For each integral topology, one can construct a, 
basis by starting vvith a parametrization of the loop momenta, of the grapls. 
Then, one needs to pick six additional propagators to get a, complete set. 
These additional propagators are important for expressing irreducible nu- 
merators. Given a, numerator, one can express this uniquely into the basis. 
The choice of numerators is by nature somevvhat arbitrary. If one focuses 
solely on additional propagators vvhich are products of differences of tvvo 
momenta, then this can be done by trial and error. A more systematic 
approach is needed, if the aim is to include graph symmetry. 

As an example, let us consider integral topology 26. Tts propagators can 
be parametrized as: 


Is, İs, ,13, (ela 1.15)”, (İs — A)”, (elə “.pi)”, (“ls pə)”, 


(“r5-16)), (Cl3Fpi-rpə)?, (13—ladl6—pə)?, (“l3--—l5—ls-pi-Epə)” . 
(7) 


This topology has tevo independent permutation symmetries, as vvill be dis- 
cussed later. It can be checked that the folloving six additional propagators 
parametrize all irreducible numerators: 


(13 — 5)”, (la — 16)”, (İs — 6)”, (a -pu”, (a — pə)”, (ls — pə)” . (8) 


In section İAppendix El choices of propagators are explicitly given for all 34 
topologies. 

A useful step tovvard computing the Feynman integrals is the so-called 
o: representation, see e.g. . The result is an integral over o: parameters, 
one for each propagator, 


1 7 : 
Ha... m) xm İ don ...don VHaş” TU-PeiF”U, (9) 
İL Tleil./0 1 


vrhere the normalization constant, vvhich is unimportant for our purposes, 
has been suppressed. U and F" are certain polynomials of the o: parameters 
of homogeneity EZ and L.-1, respectively. See m for further discussion and 

2) for a Mathematica code to compute them from given propagators. As is 
customary, the labels on the oz, parameters correspond to the position in a 
given ordered list of propagators, such as the one in equation (7) for integral 
topology 26. Equation (9) can cover negative indices, that is, numerators as 
vvell, 


2.9. Graph, symmetryi 


Several of the graphs in the set under study have a graph symmetry: 
the graphs are mapped to themselves under a, permutation of some of the 
edges and external legs. Exploiting graph symmetry can be of great help in 
simplifying and cross-checking computations if it is present. 

m İzi, graph symmetry vvas indicated briefly by the symmetry factor. 
These symmetries can be checked or found by using the built-in capabilities 
of Mathematica. These permutations can include permutations of the exter- 
nal legs. An alternative strategy is to compute the U and F” polynomials of 
the o: parameter integral, equation (6). Then, one can simply check all pos- 
sible permutations and see if they leave these polynomials invariant. Due to 
the factorial grovth of this problem, the feasibility of this analysis strongly 
depends on the algorithm: a, naive implementation ran out of steam for ten 
propagators. A faster algorithm is presented in bəl. 

A symmetry of the graph of an integral can correspond to a symmetry 
Of a, corresponding Feynman integral. For a given parametrization of the 
loop momenta of the graph, a permutation symmetry of the graph may be 


undone by a, shift of the loop momenta, combined vvith a, permutation of 
the external legs. If the latter leaves the integral invariant, the result is a 
symmetry of the Feynman integral. 

For the example of topology 26, the tvvo independent permutations are 
generated from 


ls € ls, (6: ——işebpi: 
A — —l4 İs İs, İs — —ls6 — pə, 

d 10 
İş — —İş Epi kpə, i lA — İla —la, -. 
Pı €” 2, Pi € 2, 


as can be checked explicitly in equation (7). Since the integrals, by dimen- 
sional analysis, only depend on (pı “ pə)”, they are left invariant by the 
permutation of pı and pə. 

The momentum map given above corresponds to a. simple permutation 
pattern for a given list of propagators. Let us consider, for instance, the 
scalar integral in topology 26, vrith the propagators as listed in equation 
(7). The tvvo independent permutation symmetries of the graph correspond 
to the follovring eycles: 


(UL 7), (2:8) 5,6), (5.11), (9, 120), 0) 
(0,2) 43.9) (4,10), (6,12), (7.8). 


Together vvith the trivial permutation and their product, these form a four- 
element representation of the permutation group. A given propagator can 
be mapped to other propagators in the set, but only if they are in the orbit 
of the given permutation. "he orbits of the above permutation group are 


111, 2, 7, 8), 13, 6, 9, 12), 14, 10), (5,11)). (12) 


"This can be used to simplify the sector decomposition method employed in 
FIESTA, for instance, as vvill be shovmn in section 

Given a permutation symmetry of a particular topology, ideally, addi- 
tional propagators vvould be found vhich yield a, complete basis of propa- 
gators and also incorporate the permutation symmetry. One can search for 
symmetric numerators systematically. Permutation symmetries either leave 
propagators invariant or interchange tvvo propagators. Hence, the six prop- 
agators to be added to the set contain either zero, one, tvvo or three pairs of 
interchanges and, consequently, six, four, tvvo or zero invariant propagators. 
Starting from a, generic linear combination of all four loop momenta and tvvo 
external momenta, one can construct the most general pair combination as 
vvell as the most general invariant propagator in terms of a number of free 


variables. Since propagators are only determined up to a sign by the appear- 
ance of the square, there are tvvo times tvvo possible (pairs of) polynomials 
vrhich come out of this. Hence, all possible appearances of the permutation 
symmetry can be parametrized. 

Then, one checks if the resulting set of eighteen propagators is linearly 
independent. Tn several cases (e.g. topology 20), no set exists at all, In 
addition, there are integrals vrith multiple permutation symmetries vrhich 
do not permit a choice of six additional propagators vvith all permutation 
symmetries manifest. In cases in vvhich there are sets of numerators vvith 
explicit permutation symmetries, typically many free variables exist. These 
may then be fixed by hand aiming for propagators as simple as possible. 
Care should be exercised not to choose parameters such that accidental 
Hinear dependencites in the total set of propagators are introduced. 

For the example of topology 26, the momentum maps in equation (10) 
do not map the additional numerators in equation into a, permutation of 
themselves. The folloving set of additional numerators has explicit permu- 
tation symmetry: 


(ls -- 215 — 3p)7, (İs — 215 — pə)”, 
(“iş € 216 4 pə)”, (1: — la bl: — pə)”, (13) 
(1: — 2l6 — pi — 2pə)”, (la — ls — pı pə)”. 


This example shovvs a general principle: the price of manifest permutation 
symmetry can be a, much more complicated numerator structure. İt vvould 
be very interesting to fix guiding principles for choosing sets of numerators. 

The considerations above result in tvvo sets of propagators. For each in- 
tegral topology, there is a, choice of eighteen propagators vvhich can be used 
to express amy integral in this particular topology. ln particular, any numer- 
ator can be expressed in terms of the basis. In the first set, the aim vvas to 
include numerators as simple as possible, vvhile, in the second, permutation 
symmetries vvere taken into account as much as possible. Obtaining the 
symmetric set is the result of a, fairly lengthy computation. "The simplest 


choice in the first set is listed in section (Appendix Cİ). AlI results reported 


here vvere obtained vvith this set, unless explicitly mentioned othervvise. 


3. IBP reduction and its output 


The set of Feynman integrals in equation (6) is over-complete, as there 
are relations betvveen different Feynman integrals. A particular example are 


the IBP relations İz, işl, vvrhich follov from 
D D Ö 
2— 527” ərə (14) 


VVorking out the right-hand side leads to a linear relation of different Feyn- 
man integrals. By solving a system of such equations, one may express a, 
general Feynman integral in terms of some basic integrals. 

A different vvay of phrasing the problem is to envision the system of 
IBP relations as a, large matrix equation, vvith the integrals combined into 
a, vector. "The standard vvay of solving a, problem of this type is Gaussian 
elimination. Hovvever, the output of this algorithm for a non-invertible ma- 
trix problem, such as the one under study here, depends on the ordering 
of the integrals in the vector. "Vhis is the essence of Laporta/s algorithm 
for IBP reduction : one picks an ordering of the integrals. "This or- 
dering should be such that complicated integrals are expressed in terms of 
simpler ones in general. For instance, integrals vrith smaller values of the 
parameter s (measuring numerator povvers) are preferred, as are integrals 
vrith less numbers of different propagators. 1n addition, smaller values of 
the parameter £ (measuring propagator povvers) are to be preferred. VVithin 
these general choices, in practice, the exact criteria can differ betvveen tvvo 
different implementations of Laporta/s algorithm. The results of these tvvo 
implementations are, therefore, in general not the same, but related by a, 
change of basis. 

Given a, system of equations, the output of Laporta/s algorithm is a 
reduction of all integrals in a given set dovmn to integrals vvhich cannot be 
further reduced from the given equations using the given ordering. The left- 
over integrals are knovrmn as master öntegrals. ""hese integrals depend on the 
given set of equations, although typically, if the given set of equations is 
large enough, the set of master integrals tends to converge. 

Various private and public implementations of Laporta/s algorithm exist, 
such as ATR İ55) İ55), FIRE “id El İsi and Reduze İı0), od bəl, See LiteRed İso), is1İ 
for an altermatiya approach to IBP reduction. VVe explored FIRE, Reduze and 
LiteRed in some detail for the four-loop form factor problem. Only Reduze 
succeeded in solving the problem, probably due to its better parallelisation 
implementation. 


o.1. Relevant implementation details 


Reduze is designed to run in massive parallel mode on a cluster. As an 
input, it takes a family of integrals. "Phis family is treated as an ordered 
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set: lovver integrals in the set are considered simpler. Integrals are reduced 
as far as possible to lovver integrals. "This cuts dovn on the computational 
vvorkload considerably, especially in a, large-scale problem such as the one 
under study here. As parameters, Reduze takes a range of propagator povvers 
(the £ variable) as vvell as a range of numerator povvers (the s variable) and 
constructs out of these all possible IBP relations for a, given sector. These 
are then solved sector by sector, starting from the simplest sector, vvhich is 
the bottom sector of the first member of the integral family. The number of 
relations tends to grovv very quickly along either the £ or the s axis. 

A problem vee encountered in the public version of Reduze in parallel 
mode vvas that it tended to erash vvhen handling many large files on the 
cluster file system during the identity generation stage. Due to this prob- 
lem it is unfeasible to obtain an IBP reduction vvith the public version of 
Reduze, and a fix is needed. The problem vvas that a file vvas reported as 
not present in the file system, vrhile a, manual check thereafter did uncover 
this particular file. "This is due to the internal file handling structure of 
Reduzel” vrhere all processes are allovved to read from and vrrite to disc. As 
a vvork-around, the program can be forcedİ to sleep for a, time-out of fifty 
seconds, follovved by a, 500-second time-out if the file is still not found. This 
tvvo-stage timeout resolved the crashing problem due to disk space usage. 
The loss of productivity for a fevv minutes is a small price to pay to prevent 
the total collapse of the computation. 

During the running of the four-loop reduction, Reduze still crashed oc- 
casionally vrhen handling integrals in sectors vvith many, typically tvvelve 
or eleven propagators. These integrals require a large memory usage each, 
vrhich can easily overvvhelm one of the nodes of the cluster. Sometimes this 
can be due to bad scheduling of processes, vvith all zoof processes vvhich con- 
sume most memory running on a single node: this situation can somevhat 
be avoided at start-up by passing scheduling instructions to the message- 
passing-interface (MP1) program? By contrast, for lovver propagator count, 
the limiting factor for speed is mostly CPU time. These tvvo situations re- 
quire opposite numbers of allovved parallel processes, vvhich in turn requires 
occasional input by the user during running. It should be noted that, al- 
though Reduze vvorks in parallel, there is a, limit to the number of cores 
vrhich are assigned to a, single process vvhile still increasing computational 


SVVe thank A. von Manteuffel for explaining this. 
“This can be done in the function set-)ob.-status in the file Pobqueue.cpp. 
5One uses -bynode in openmpi. 
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speed. This saturation is sector dependent, but certainly under a, hundred 
COT€S. 

The chosen İBP relations are in the range £ X€ 12 and Ü € s € 2, vith 
some extensions to £ — 13 for specific integral topologies. Since Reduze 
features a, choice of used relations, there can be problems vvith the reduction 
at the outer edges of these relations: there can be unresolved integrals, A 
simple check of this is to observe the file size of the reductions in a sector 
and compare it to the average for other integral topologies in a, sector vvith 
the same nuümbers of propagators. A sudden increase in file size tends to 
indicate unreduced integrals: the unreduced integrals appear vvith much 
more complicated prefactors than the others, typically involving polynomial 
ratios vvith polynomial orders an order of magnitude above the norm. This 
is more than fust a nuisance. The abnormal prefactors tend to be highly 
divergent in the limit € — Ü, typically oc ziz or vvorse. Since the form factor 
itself is expected to diverge at vvorst as na this implies massive cancellations 
among the integrals. This particular problem occurred in the basis of section 
Appendix Clin integral topologies 20 and 26, for integrals vvith at most ten 
propagators. This problem vvas resolved using consistency conditions, see 
tvvo paragraphs dovn. 

VVe reduced using the set of numerators in section İAppendix Cİ A sec- 
ondary computation using numerators vrith explicit permutation symmetries 
failed: at a, certain eleven-propagator integral, the MPI protocol ran into an 
issue vvith maximal size of a, single message, vvhich is set by the maximum 
size of an integer on a, 32-bit machine, see for an explanation of the 
problem and a possible resolution. "This is a result of a choice of protocol 
and not related to the physical architecture used. In principle, this problem 
could be circumvented by revvriting the Reduze code to handle messages of 
this size. 

Since the set of numerators vvith the simplest structures vvas used, inte- 
gral relations from permutation symmetries are not explicit. "This allovved 
a simple check on the results: take a, graph symmetry relation and use the 
IBP reduction. If Reduze found the relation, the result must be zero. VVe 
checked permutation symmetry on maximal-propagator integrals vvith up- 
to-quadratic numerators. Surprisingly, one single additional relation vvas 
found vhich involved exactly all master integrals vrith abnormal prefactors. 
Plugging this relation into the reduction eliminated the abnormal prefactors, 
reducing them to much more normal-looking ones. 
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9.2. Results 


The computation led to a result after several months of computation on 
a, fairly large cluster İİİ The complete list of the master integrals obtained 
vvith this reduction is attached, both in the planar and non-planar sectors, as 
vvell as a, sample result for topology 25. Overall, the integrals vrhich appear 
after reduction are simpler. "There are tvvo integrals vvith tvvo numerator 
povvers. One of these contains only eight propagators and is easy to integrate 
numerically. "The other, hovvever, has the full tvvelve propagators and occurs 
in integral topology 26. It turns out that this particular integral can be 
integrated numerically (see section /5.1). 

Any basis of master integrals listed is arbitrary to some extent. For a 
given sector, an obvfous first choice for a master integral is its corner integral. 
Beyond this, there is a choice of a, (typically single) numerator or doubled- 
up propagator. The statistics in Table Hİ relates to the direct output of 
Reduze, taking into account the one extra, relation found from permutation 
symmetry. The statistics is split into those integrals vrhich contribute to the 
planar and non-planar parts of the form factor. 


Table 1: Master integral statistics of obtained IBP reduction. 


(a) planar form factor (b) non-planar form factor 
sl s—2 sl s—2 
6 0 10 1 
2 0 S) 0 
9 0 10 0 
1 0 1 0 
4 1 3 1 
0 0 0 0 
0 0 0 0 
0 0 0 0 
22 1 27 2 


In these tables, there is quite some double counting of integrals vvhich 
appear both in the planar and non-planar sectors. Using the knovn result 
for the sum of the integrals in the planar sector, one can eliminate several of 
the integrals in the non-planar sector. "This is more eflicient than might be 


“Memory usage on single sectors can exceed 200 GB. More detailed information is 
available on request. 
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expected, as there are four integral topologies (21, 25, 30 and 31) vrhich each 
appear vrith the same overall color factor, 2/V5 -- 24N2. This can be used 
to eliminate these integral topologies from the non-planar sector in favor of 
the integrals in the purely planar sector and the knovn planar-form-factor 
result. The resulting sum contains more integrals (260 in total) than the 
non-planar sector by itself, 

The choice of master integrals in the provided results is taken from the 
output of Reduze. İt is possible to choose other sets of master integrals, in 
particular the ones vvith numerators, to aim at simplifying, for instance, the 
prefactors of the integrals in the form factor or to make use of a permutation 
symmetry of the scalar skeleton of the integral. In particular, one aim could 
be to make these as vveakly divergent as possible in the limit € — 0. The 
povver of divergence in this limit determines the order to vvhich the integral 
is to be expanded to get the cusp anomalous dimension. 


Cross-check: Multiple reductions 

The IBP reduction vvas run tvvice, vvyith a, second reduction involving 
a, reshuffling of the integral topologies to put planar topologies first. İt 
vvas checked that the tevo reductions commute: reducing vvith the second 
reduction follovved by the first gives the same result as reducing vith the 
first, after taking into account the extra relation found from permutation 
symmetries. Before taking into account this relation, there are seeming 
mismatches betvveen several of the integrals vvhich appear in the full four- 
loop form factor: these mismatches are all proportional to the extra relation. 
"This is a consistency check for the IBP reduction. 


Ə.Ə. Cancellations in the form factor after TBP reduction 


The TİBP reduction displays tvvo separate patterns of interesting cancella- 
tions for the form factor result: one pattern takes place vvithin the separate 
integral topologies, vvrhile the second is betvveen different integral topologies. 

The first pattern is that, for integral topologies 18 and 20, the contribu- 
tions to the form factor vanish after IBP reduction. Generically, the Feyn- 
man integrals in these topologies do not vanish: there are tvvelve-propagator 
master integrals. Hovvever, the specific combination of the integrals in the 
AN — 4 SYM form factor reduce to zero. Considering the complicated nu- 
merators for topologies 18 and 20, this is rather nontrivial and not obvious 
at all at the integrand level, It should be noted that these integrals already 
do not contribute to the form factor: their color factors are also zero. Such 
cancellations, hovvever, do not happen for the other four zero-color-factor 
topologies, 8, 11, 15 and 16, vvhose integrals reduce to nontrivial results. 
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It vvould be very interesting to find an explanation for these vanishing in- 
tegrals, as, in physics, a zero almost alvvays has a, physical and important 
symmetry explanation. 

Another intriguing pattern emerges vvhen tracking the dependence on 
the single free parameter in the form factor integrand vvhich vvas left in 
Hs). There, the free parameter in the form factor multiplies a, combination 
Of 24 out of 34 different integral topologies vvhich obeys all the physical 
constraints imposed in separately. 1n particular, this combination is 
color-kinematic dual by itself, One may expect that, by applying other, 
more complicated cuts vvhich are not exploited in Hs), this parameter may 
be fixed uniquely. The surprise is that, vvhen summed vrith the correct color 
and symmetry factors (and only then), the dependence on the free parameter 
drops out of the form factor after IBP reduction. This means that the color- 
kinematic-dual representation of the four-loop form factor contains a truly 
free parameterl 
.— This detailed cancellation involves an order of fifty different master in- 
” R: tegrals vvith up to ten propagators and occurs in both the planar and non- 
.. planar sectors. Most of the cancellations of the coefficients of the master 
integrals involve integrals out of all of the eighteen integral topologies vvhich 
involve the free parameter and have nontrivial color factors. As an intricate 
example, consider the coefficient of the master integral 


INTİ2, 1-1, 0, 1,1, 1, 1,1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0)1, (15) 


vrhich occurs in topologies 12, 14, 15, 17 and 19. The coeflicient in topology 
12 reads: 


.— —242 H 2861e — 4906€? — 37519€” -- 133706e? — 18520c? — 378200€6 -- 366000e7 
F” 2164(1 H 4e)(—1 4 5e)(1 4 5e) 
(16) 
vrhich is exactly —2 times the coefficient of this master integral in topology 
19. Taking into account the symmetry factors of topologies 12 and 19, 
vrhich are 2 and 1, respectively, this leads to a, cancellation. Similarly, the 
coeflicients of this master integral in topology 14 cancels against that in 
topology 17, vrhile the color factor of topology 15 is zero. The other master 
integrals typically involve much more intricate cancellations. 

It vvould be very interesting to reach a, better understanding of the un- 
covered cancellations. On a practical level, this pattern completes the de- 
termination of the form factor integrand İl. To our knovvledge, this is the 
first example of a, color-kinematic-dual integrand representation that con- 
tains a, free parameter. It vvould be fascinating to see if there is any deeper 
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physical meaning of the free parameter, for example, as due to some gauge 
symmetries. This free parameter can be set to any value at the start of the 
computation. A natural value is —1, vvhich simplifies the naive appearance 
Of the integrals. "This cancellation provides strong evidence for consistency 
of both the integrand result in and the TBP reduction performed in the 
present paper. 


4. Counting master integrals using algebraic techniques 


m İo1l, Lee and Pomeransky developed a, method to count the number 
of master integrals using a technique adapted from algebraic geometry. "Phis 
is partially based on previous vvork by Baikov İs3). This method only relies 
on the topology of the integral and, therefore, applies for arbitrary values 
Of £ and s. İn particular, it is independent of obtaining an explicit IBP 
reduction and, hence, of the results of the previous section. This idea vvas 
also implemented in a Mathematica package called Mint . İn this section, 
vve develop and apply this method to count the number of master integrals 
of four-loop form factors. 

Similarly to the start of TBP reduction, one first needs to determine the 
unique physical sectors for a, given integral topology. As described in İ:l, 
it is convenient to perform this step using LiteRed İso, İsi. mm particular, 
LiteRed determines the sector symmetries, vrhich can be taken into account 
automatically in Mint. İn the next step, the counting of master integrals is 
done independently for each physical sector. 

Let us consider a physical sector of £-loop topology and r propagators. 
From the associated polynomials U and F” in equation (9), one can define a 
nevv polynomial?l 

G(a) — U(d) — F(a), (17) 


vvrhich has polynomial degree AL -F 1. Using the duality of homology and 
cohomology (see for more details), the number of master integrals is 
related to the number of independent integral contours. The latter depends 
on the analytic structure of the polynomial G and is characterized by the 
so-called proper critical points. "These are the stationary points at vvhich G 


TEquivalently, one can consider a different polynomial using Baikov”s representation 
as discussed in in ss). For all cases being checked, the tvvo polynomials alvrays give 
the same result, vvhile the use of the polynomial G is usually much more efficient in the 
computation. 
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does not vanish and are, thus, defined as the solutions of 


əs 
Öo:, m" 


1575 and G z£ 0. (18) 


Proper critical points may or may not exist for a given G polynomial, If 
there is no critical point, this implies that the corner integral of the sector 
is reducible, so that this sector does not contribute a, master integral. On 
the other hand, if critical points do exist, the corner integral is irreducible 
vrith respect to IBP. This condition of irreducibility vvas first shovmn in : 

A proper critical point may have a malziplicity vvhich is measured by a, 
topological invariant, the so-called Milnor number at the critical point. "The 
number of master integrals is the sum over the Milnor numbers /M, of all 
proper critical points, 


zi of master integrals — x M:, (19) 


?€proper critical points 


provided that the proper critical points are isolated. VVe vvill explain the 
meaning of ?solated eritical point shortly. Im the case of non-isolated eritical 
points, more vvork is required , See also an explicit example in section 
belovv. 

Thus, the problem of counting master integrals is reduced to the one 
of computing the critical points of a system of nonlinear polynomial equa- 
tions given by equation (18), vrhich is usually highly nontrivial to solve. 
Fortunately, one can apply povverful algebraic approaches to make this ?ob 
simple, based on the so-called Gröbner basis technique. Readers vvrho are 
not familiar vrith the Gröbner basis and related concepts may consult 
for a, pedagogical introduction. Here, only the basic procedure is outlined. 

"The solution space of a system of polynomial equations is called affine 
variety associated vvith an ideal determined by the equations. For the system 
of equations (18), the ideal can be defined as 


olel öG 
1—0——..:— —1 20 
(əə ” Om , ooG . , ( ) 


vrhere, in the last term, one introduces an additional parameter op vvhich 
forces the polynomial G to be non-vanishing at the eritical point. "The 
solutions of the affine variety can be obtained by computing the Gröbner 
basis of equation (20), 


gb(1) — 19g1,92:...:9kl . (21) 
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Once the Gröbner basis is obtained, it becomes relatively trivial to find the 
solutions. In our problem, vve only concern ourselves vvith the number of 
solutions. "This, in practice, can be conveniently obtained by counting the 
number of irreducible monomials in the Gröbner basisl 

The above-described procedure is implemented in Mint and applies 
straightforvvardly to many simple cases: for example, applying it to the 
three-loop form factor, vve obtain the result summarized in Table The 
numbers of master integrals are consistent vvith the reduction of LiteRed. 


Table 2: Master integral counting for the six three-loop form factor topolo- 
gies given in Figure 5 of İH). 


topology (1) (2) (8) (4) (5) (6) 
gi Mis via Mint İ 9 10 10 10 14 10 


For the four-loop form factor, hovvever, one encounters tvvo problems 
vrhen trying to apply the Mint package. The first problem is that the set of 
critical points can form an affine variety of dimension S 1. "V"his corresponds 
to the non-isolated critical-point case. Mint cannot deal vvith such cases 
automatically in its present version, but denotes them as /ndeterminate. 
The number of such cases for the integral topologies at hand are summarized 
in TableBl Although the absolute number of these cases is lovv, they tend to 
occur in the more complicated sectors. "VPhey require some further vvork as 
described in . VVe employed a similar procedure and present an explicit 
example in subsection İ4.1l 


Table 3: /Tndeterminate sectors from Mint for the 34 topologies given in 
Figures 8 and 9 of l3). 


topology (14) (15) (16) (17) (18) (20) (22) (23) (24) 
z£ mdeterminate 1 1 1 1 1 2 1 1 1 
topology (25) (26) (27) (28) (29) (30) (32) (33) (34) 


z£ mdeterminate 1 3 1 1 3 1 1 1 1 


Another problem is that, in quite a lot of cases, Mint gets stuck at the 
step of Constnucting Gröbner basis .... "This is due to the complexity of four- 


S More precisely, one needs to compute the reduced Gröbner basis, vvhich is unique for 
a given monomial ordering, see e.g. 64). The number of irreducible monomials also takes 
the multiplicity of eritical points (1.e. the Milnor number) into account. 
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loop integrals for vvhich Mathematica cannot finish the computation of the 
corresponding Gröbner basis. There are many knovrn other packages vvrhich 
can compute a Gröbner basis more efliciently. VVe applied Macaulay2 İs), 
vrhich can be conveniently used in Mathematica through mathematicaM2.m 

, and vve employed Singular in several most complicated cases. 

Based on Mint and vith some extra, effort to solve the above tvvo prob- 
lems, vve obtained the numbers of master integrals for the 34 different topolo- 
gies as summarized in Tablel4l 


Table 4: Master integral statistics: Mint results for the 34 topologies given 
in Figures 8 and 9 of İH). 


topology İ (1) (2) (8) (4) (5) (6) (7) (8) (9) (10) (11) (12) 
q MIs 28 35 3 34 ə9 48 T1 52 55 69 52 95 
topology 1 (13) (14) (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) 
qF MIs 76 öT 95 93 106 g4 105 39 45 66 41 T5 
topology İ (25) (26) (27) (28) (29) (30) (31) (32) (33) (34) 
qF MIs 55 Ts 92 69 93 94 Əs ə əl ə9 


Given the number of master integrals, Mint can also generate a set of 
master integrals. ln most cases, a sector contains only one master integral, 
in vvhich case the corner integral is chosen as the master integral, If there 
are more than one master integrals in a sector, Mint suggests integrals vvith 
higher povvers of propagators (vvhile Reduze typically chooses integrals vvith 
numerators). Hovvever, the choice of a double-propagator integral is rather 
heuristic: therefore, to find out vvhether it is a genuine master integral re- 
quires further checks. In our case, vve could check them by using Reduze, as 
vvill be discussed shortly. 

"To obtain the master integrals for the four-loop form factor, vve still need 
to combine the master integrals of Tablelltogether. Obviously, same master 
integrals may appear in different topologies of Tablelil This may be detected 
by comparing their graph topology: if tvvo graphs are isomorphic to each 
other, they are the same integral. After this step, vve find that there are 
only 280 non-isomorphic master integrals: 244 vvith only simple propagators 
plus 36 containing a double propagator 


9Note that vve combined all 34 topologies including six zero-color-factor topologies: (8), 
(11), (15), (16), (18) and (20). If vve exclude these six topologies, the number is reduced 
to 267, vvhich is enough for the /V/ — 4 SYM form factor. But for QCD vvith fundamental 
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Here, one should be cautious that there could be relations betvveen inte- 
grals vvhich are apparently not due to graph isomorphism, see e.g. İz0, for 
some explicit examples. Such relations can be found by considering Feynman 
parameter representations of integrals and checking if the tvvo representa- 
tions are related to each other by a permutation of Feynman parameters. A 
fast algorithm for such a check vvas proposed in . Alternatively, one can 
also consider the so-called grap) matroid technique to detect such hidden 
relations betvveen different integrals, as implemented in Reduze . For us, 
there is an easier vvay to detect such relations, namely applying Reduze to 
reduce the basis integrals obtained by Mint. VVe checked that all 244 corner 
integrals correspond to basis integrals as determined by Reduze. Thus, they 
are indeed master integrals. 

As mentioned above, the choice of master integrals vvith double propaga- 
tors by Mint requires further checks. VVe applied Reduze to the 32 double- 
propagators integrals vvhich do not appear in the tvvelve-propagator sector of 
topology 26 and found that four of the suggested master integrals are actu- 
ally reduciblellil "This does not change the number of Mint master integrals, 
but simply means that a different integral from this sector must be chosen 
as a, master integral. By changing the position of the double propagator, vve 
indeed found irreducible integrals. The tevelve-propagator sector of integral 
topology 26 is analyzed further belovr. Thus, vve obtained the final list of 
280 genuine master integrals. 

VVe vvould like to emphasize that the counting of master integrals based 
on Mint only relies on the topologies of given integrals and applies to ar- 
bitrary numerators, therefore, it applies to any theory, including QCD. VVe 
remind that the counting of Reduze given in Table İİ concerns particular 
master integrals appearing in the reduction of the /V — 4 SYM form factor 
result. 


An interesting mismatch, betuveen Reduze and Mint 


Some caution is needed vvhen assessing the output of Mint: this package 
is knovn to under-report the number of master integrals in very rare cases. 
To check this, vve collected all master integrals appearing in the IBP 
reduction by Reduze, vith £ — 12 and s — 2 and vrith £ — 13 and s — 1. It 
turns out that Reduze yields more master integrals. Many of them consist 


quarks, one likely needs all 280 master integrals. 

1ÜNote that Reduze tends to choose integrals vvith numerators rather than double prop- 
agators. So, örreducible means that, in the reduction result of Reduze, there is an integral 
vvith the same propagators containing a numerator. 
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Of tvvelve propagators and one to several numerators. They are at the edge 
of the reduction setting, and it is very likely that they may be reduced by 
including more IBP identities. "These vvill not be discussed further here, 
details are available on request. VVhat is unexpected is that there are nine 
corner integrals vvith only ten or less propagators, for example 


INTI2 10,0.1,1.1,1,1:1.0.0,0-1.0.0.1:0,0,0h: (22) 


vrith only eight propagators. Mint shovvs that they contain no proper critical 
points, vvhile Reduze and FIREB take them as master integrals. Due to the 
lov number of propagators, most related İBP identities should already be 
taken into account by Reduze. Therefore, it is likely that these integrals 
are truly master integralsİLİİ "This vvould mean that the counting method of 
Mint requires some modification, althougl it is correct in the vast ma)ority 
of cases. VVe cannot rule out, hovvever, that there are integral relations 
vrhich render this set of nine corner integrals reducible. It vvould be very 
interesting to understand this issue and see vvhich is right. 

m addition, there is one more master integral in the Reduze basis, vvhich 
appears to be an artifact of the vvay vve solved the permutation symmetry 
consistency condition. 


4.1. Tuvo urorked-out eramples 

Here, vve give further details by considering tvvo examples vvhich cannot 
be computed directly by Mint. 

Let us first consider an example corresponding to an indeterminate case: 
a subsector of topologiy 26 vvith eleven propagators, 


INTİ26, (1,1,1,0,1,1,1,1,1,1,1,1,0,0,0,0,0,01. (23) 


From G £— U H.F, vhich are polynomials of oq (4 — 1,...,11) and of ho- 
mogeneity five, one can construct the ideal 7 defined in equation (20). The 
corresponding Gröbner basis cannot be computed by Mathematica, but can 
be calculated using Macaulay2 . “The quotient space of the resulting 
ideal Z — GB(T) turns out to be infinite dimensionalİ"”İ One then computes 
the primary decomposition of 7 (e.g. using Macaulay2), vrhich gives 


T—nnTə, (24) 


T.Some other examples of tvvo-loop integrals shovring mismatch betvveen Mint and IBP 
reduction vvere also found by V. Smirnov İs. 

12The dimensionality of the quotient space is given by the number of irreducible mono- 
mials, vvhich can be directly counted from the reduced Gröbner basis as implemented in 
Mint İ21). 
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vvhere 


T) —(186624op — 3125, aq, 5az “12, 5oş — 12, 5aq —- 12, 5aş — 12, o, Boy -- 12, 
Bag “b 12, 5ag — 12, 5oqo -F 12, 5aqi 12), 
Tə — 7290p — 16, 20 “EF 20 —£ 3, az 20) FE 6, aş — 3, 2044 — 20ag “3, üş — 3, 
206 — 20aqi 1.3, öz 20“ 6, 20) aş F 30 EE 30 “- 9) 209 — 3, 2040 — 2ai 3). 


The quotient space of the first ideal is one-dimensional, vvhich is obvious, 
since there is only one solution of the corresponding polynomial system. On 
the other hand, the quotient space of the second ideal is infinite dimensional. 
This can be treated as the example discussed in İl. Rather than consider- 
ing the full ideal 7ə, one can simply consider the variety determined by the 
equation Ğ(as, on?) — Ü, vrhere 


G(as, 6n)) € 20niaiş FE 30 EF 3oni 9 (25) 


is the only quadratic element in 7ə. "TThen, one computes the dimensionality 
Of the quotient space of the ideal 


İ — (9Ğ Əos , OĞ /Öoqi : oo — 1) : (26) 


vrhich turns out to be one. Altogether, one finds that the number of bases 
is tvvo. 

As another illustrative example, vve consider the sector of topology 26 
vvrith full tvvelve propagators. In this case, the computation of the Gröbner 
basis turns out to be hard even for Macaulay2. In this case, ve apply 
Singular, vrhich solves the problem in less than a minute using the method 
slimgb. "The quotient space of the resulting Gröbner basis 7 — GB(T) is 
seven-dimensional. One still needs to consider the symmetry vrhich may 
reduce the number of independent master integrals. 1n order to make the 
symmetry obvious, one can compute the primary decomposition of Z, vvhich 
gives 


TZ — Tı nT7ə 113 174 17: 125, (27) 
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vvhere 


1) — 729ap — 16, 4o HE 3, az — 8, oxa — 3, iq, 40iş, -E 9, dış — 3, oy — 3, 20iş HE 9, 4ag 3, 
20x0 — 3, oq “6, 202 FF 9), 

Tə —7290p — 16, o — 3, 20 F 9, o — 3, diq, Ox -E 6, dig — 3, 4oyy —E 3, aş — 3, 2a 9, 
20x0 — 3, 40 F 9, 402 3), 

İş —(729ap — 16, 20 € 9, az — 3, 20iş 9, 204 — 8, ox -t 6, 4ig F 3, öz — 3, 4oiş 3, 
ag — 3, aqo, 4ani 9, aqə — 3), 

T4 — 72900 — 16, oz — 3, 4aə F 3, 4oxş -E 3, 204 — 3, 40) -E 9, 2aig “9, 20? -E 9), dış — 3, 
og — 3, Oq0) Oq 6, ış — 3), 

Ts —256ag “1, oq “2, az 2, aş 2, 04 — 2, aş 6, aç “2, öz -E. 2, dış 2, ag 2, 
ön0 — 2, önü — 6, qə “.2), 

T6 —(65536co0 — 3125, 5aq “- 5aqz — 8, 5ozə -E 5orə — 8, dış — o?) 20c4  5a2 — 16, 
Əqiş -b 16, aş — qə, Süz -E 5a — 8, 5aiş — 502 — 8, diş — də, 
200xio “EF 5oriə — 16, 5or “EF 16, 25o2,) — 400 — 128). 


ı. The quotlent space of the last ideal is tvvo-dimensional, vvhile the others are 
— all one-dimensional. So, naive counting vvould give seven master integrals. 
-. Topology 26, hovvever, has permutation symmetry, see equation (LI). Tt is 
not hard to see that the first four ideals are all related to each other by this 
symmetry:, therefore, they contribute only one independent master integral. 
m total, one finds that there are four independent master integrals in this 
sector. For this integral topology, Reduze picks four tvvelve-propagator inte- 
grals as master integrals: the corner integral, tvvo linear-numerator integrals 
and one vvith a quadratic numerator. Hence, the number of master integrals 
matches betvveen Mint and Reduze. 


5. Tovvards numerical integration 


Given a, set of master integrals, the next step is to perform the inte- 
gration. A significant part of the master integrals obtained in the previous 
sections are very complicated and seem to be on or slightly over the very 
edge of vvhat current integration methods can achieve. There are various 
methods to perform numerical integration in the literature, see for an 
OVETVİEVV. 

It should be said that, in principle, Mellin-Barnes (MB) representations, 
vrhich can be used very efficiently in solving planar topology integrals 169, İd, 
İr1), , may be the optimal vveapon of chotce for four-loop integrals, if the 
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difficulties in obtaining efficient MB representations for non-planar integrals 
vrithout hidden singularities can be resolved. See İ9), page 124, for a brief 
description of the problem. Progress in obtaining up to and including some 
three-loop non-planar representations in an automatic vvay vvas reported in 
Vəl. This may be used, at least in principle, to obtain MB representations 
for those integral topologies containing either a, bubble or a, triangle vvith 
tyvo on-shell legs. 

VVe mostly explored numerical integration using sector decomposition 

. ln general, this is a method to turn a highly divergent integral into a 
sum over finite integrals, by resolving the singularities of the Feynman pa- 
rameter integral, "These finite integrals can then be integrated, for instance, 
using the CUBA library Vəl. VVe are avvare of three different implementations 
of this idea in public codes: sector.decomposition Val, SecDec İzA İrs, İrə) 
and FIESTA İsd, İsıl, isə). Only the latter program in its latest versions (3.x) 
vvas able to handle some of the most difficult integrals of the problem at 
hand. 


ö.1. Relevant implementation details 


There are various strategies for resolving the singularities of a, given 
Feynman integral. FTESTA offers a choice of these. Most relevant here are 
strategies S and X ral. Strategy KU İs3), vrhich, in principle, is highly 
efficient in terms of sectors, did not resolve any of the sectors of the most 
complicated integrals (the corner integral of topology 25) in either of its 
three variants and vvas not pursued further. Strategy X, if it terminates, 
tends to yield better results in terms of the number of sectors than strategy 
S. 

FIESTA splits the Feynman parameter integral in equation (9) into dif- 
ferent pr?maryı sectors. 1n each primary sector, one particular Feynman 
parameter is smaller than the others. If the integrand is invariant under a, 
permutation of the o: parameters, then these primary sectors may be related. 
A particular permutation symmetry of the o: parameters descends from the 
parent graph. This symmetry may then be used to reduce the computa- 
tional vvorkload. Let us recall the symmetry property of the scalar integral 
of topology 26 given in equation (11). From its orbits of the permutation 
group in equation (12), 


1112, 7,81 13,6,9, 121 14, 10), 15,117). (28) 


vve can see that, instead of computing primary sectors 1, 2, 7 and 8 sep- 
arately, one better yust computes one of these, as the others vvill give the 
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same numerical ansvver. In practice, it can matter vvhich coefficient is picked, 
since the structure of the integrand for the o: parameters is different: this 
can result in different numbers of sectors. 

The integration time varies: for integrals vvith tvvelve propagators, the 
preparation time alone is typically measured in days to vveeks. By contrast, 
all integrals vvith less than eleven propagators are highly tractable using 
FIESTA, typically requiring less than a day. 


Quadratic-mumerator integral in topology 26 


The quadratic-numerator integral in topology 26 is naively expected to 
be one of the hardest integrals in the full set. It can be integrated using 
FIESTA. Integral topology 26 has tvvo permutation symmetries, see equation 
(10). To make use of this, one needs to pick a symmetric quadratic nu- 
merator. "That is, one needs a, numerator vvhich is mapped to itself under 
both symmetries and vvhich is related to the quadratic-numerator master 
integral in the output of Reduze by IBP reduction. Moreover, to facilitate 
the computation of the double derivative in the o: parametrization, it is ad- 
vantageous to pick a numerator that is linear in the loop momenta., İt is 
“ə easy to check that 
— (ls - (bi — pə)” (29) 
transforms into itself under the symmetries in equation (10). VVe checked 
that this is related to the quadratic master integral from the explicit IBP 
reduction. VVith the symmetries in place, one only needs to compute four out 
of tvvelve primary-sector coeflicients in FIESTA. VVe picked the coefficients 
14, 0, 4, 2, 2, 0, 0, 0, 0, 0, 0,0) for the order of propagators in equation fZ), 
. — basically choosing the lovvest from the orbits in equation (28). 

The resulting integral can be integrated vvith some numerical effort, giv- 
ing the result 


INTƏel(ls - (pı — pə))”) —(—0.032986 -E 2. 1077) €” -£ (0.0694456 z- 1.01 - 107”) €”7 
-- (1.3506 -- 0.0002) €”” —E (—2.68804 -k 0.00317) €”? 

-- (-6.23707 -£ 0.04013) €”?  (12.6763 -k 2.0782) €” 

-- (1234.49 -£ 32.97) €”? 1 O(e€”1). (30) 


Preparing for the numerical integration takes about five days. "The inte- 
gration time itself strongly depends on the integration parameters in the 
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Vegas integration algorithm of the CUBA libraryll“l VVe checked by running 
vvith different numbers of evaluation points that the first five coeflicients are 
numerically stable under variations of the relevant parameters, especially 
maxeval. In other vvords, it is likely that their errors are in the reported 
ranges. 


6. Summary and discussion 


ln this paper, vve obtained an explicit integral reduction for the four- 
loop Sudakov form factor in /V/ — 4 SYM based on the integrand result in 
İLİ and by explotting the Reduze package İv). This reduces the challenging 
four-loop problem, in particular the computation of the unknovn non-planar 
cusp anomalous dimension, dovn to the evaluation of an explicit set of given 
master integrals. 

VVe also obtained master integrals by using the algebraic method intro- 
duced in b:l. This provides a first non-trivial implementation of the method 
to a, complete observable. Limits of the Mint package and their resolutions 
are explained. "The resulting integrals are, vvith some qualifications, the 
master integrals of quite generic form factors. For instance, they are ex- 
pected to cover all integrals vrhich appear in the QCD result. By comparing 
vrith the results of IBP reduction, vve found interesting mismatches, namely, 
existing masters integrals from IBP reduction are taken to be reducible by 
the Mint package. İt vvould be important to understand the reason of this 
discrepancy. 

Based on the reduction result, vve furthermore observed surprising can- 
cellations in the form factor. First, the form factor integrand obtained in 
contains a, free parameter. VVhile it vvas unclear vrhether this parame- 
ter could be uniquely fixed by further physical constraints, vve shovved that 
the terms depending on this parameter all magically cancel upon using IBP 
identities. Therefore, this parameter is truly free and gives the first example 
Of a, color-kinematic-dual integrand vvith such a free parameter. PFurther- 
more, tvvo integrals, corresponding to topologies 18 and 20, turn out to be 
exactly zero, vrhich is not obvitous at all at the integrand level, Tt vvould be 
very interesting to have a physical understanding of these intriguing cancel- 
lations. 

VVhat remains to be done is to compute the master integrals. Hovvever, 
for many of them, this is currently yust out of reach for publicly available 


13A brief exploration of the other algorithms in the CUBA library quickly shovvs that 
Vegas is the optimal choice here, as for both time and accuracy. 
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(implementations of) integration methods. "There is an important analogue 
to the three-loop case, vvhich vve vvould like to dravr attention to. In that 
case, a, reduction to master integrals vvas obtained in lil. Their integration 
vvas first performed almost three years later in 13), using partly numerical 
methods. An analytic result vvas then reported in 84 , vvith an important 
cross-check in İ85). Currently, the three-loop form factor integrals often ap- 
pear as an example application of nevr integration and organization methods, 
e.g. in İS6, 

A nevv feature in the four-loop case is the appearance of several master 
integrals beyond corner integrals. These can be taken to be either doubled- 
up propagators or single numerators. "The latter can, in principle, be traded 
for higher-dimensional scalar integrals using the Feynman parametrization 
in equation (9). The problem is that these scalar integrals feature four dofs: 
there is a sum over integrals vvith four doubled-up propagators. VVithout 
an IBP reduction of these integrals dovn to the basis, this is an unpractical 
result. An explicit IBP reduction vvith up to four dots vvould enable the use 
Of several methods. For instance, one could use dimensional recurrences 38 
to either compute the integrals numerically or to express the form factor 
in terms of a quasi-finite basis İs7). An explicit quasi-finite integral reduction 
for the three-loop form factor has been announced 0). The needed IBP 
reductions at four loops seem very hard to obtain vvithin the current version 
of Reduze and simply out of reach for anything else on the public market. See 
for a, proposal of an improved algorithm vvhich is particularly suitable 
for single-scale reduction problems. 

There is a practical issue vvhen calculating the form factor vvhich concerns 
the choice of master integrals. "Phere are several choices of master integrals 
one could make as for both linear numerators and doubled propagators. 
VVhich one to choose should then be guided by öntegrabəölity, 1.e. by vvhich 
of these choices they can actually be integrated. Moreover, as for numerical 
evaluation, there are issues of precision and speed to be considered. One 
vvould expect doubled-up propagators to give simpler Feynman parameter 
integrals, vvhile linear numerators tend to give less divergent coeffcients. A 
canonical example of the tvvo choices are the master integrals of the double- 
box-integral topology. One can either choose a. doubled-propagator basis 
element, as first done in l92). Alternatively, one can use a linear-numerator 
integral, as in İ93). The latter has a less divergent prefactor as a, function 
of €. VVhich choice is the more practical must be decided on a case-by-case 
basis though: vve observed cases vvhere no difference in expansion order is 
present. 

"There are various other avenues to be explored vvhich may lead to a 
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full computation of the cusp anomalous dimension through the form factor 
route. One vvould be to explore MB representations of four-loop non-planar 
integrals. An obstacle here is that naive implementations of this idea quickly 
lead to high-order MB representations containing an order of one hundred 
MB integrations. Any algorithm vvhich could curb this number vrould offer a 
vvay forvvard. VVhile numerical evaluation may have the precision issue, ana- 
İytic methods vvould be more desirable. Linearly reducible integrals may be 
evaluated analytically by Hypınt g41. This vvorks particularly efficiently for 
integrals of less singular behavior 37İ. An alternative route for computing 
the cusp anomalous dimension is to extend integration techniques based on 
differential equations İ95) to form factors by making one extra. leg off-shell 
and then taking the massless limit. Controlling this limit after integration 
might also be an obstacle. Apart from that, IBP reduction, a necessary 
step in constructing the differential equation, vvould be more challenging for 
deformed integrals. See 86İ for an application to the three-loop case. 
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Appendix A. Three-loop form factor check 


The three-loop form factor vvas first computed in il. vve started from 
the expression obtained in , vvhich can be shovn to be equivalent. The 
problem consists of six integrals vvith at most a, single numerator. "The 
integrals contain nine propagators, and three additional numerators have to 
be specified to obtain a definite basis. Five of these integrals have a single 
permutation symmetry. 

IBP reduction using Reduze on a part of our local computing cluster 
running in parallel mode took about 1.5 days. The reduction yielded a set 
of master integrals, the most complicated of vvhich indeed correspond to 
those obtained in i4İ. "These integrals vvere then computed using FIESTA, 
running for a fev days vvrith their standard precision setting. The result for 
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the form factor 18 


FÖ) - 0.166606 0.00012 0.00054 1.10182 0.93788 O 1 
.— 005 DU: €i ü - ü e? .. 
(A.1) 
vrhich is to be compared to the knovrn result, e.g. from bl, 
1  11((3) 247v3 1 
m "o. -Oİ-). A.2 
— nurlulüur”r.uub”” a .. 


The difference betvveen the exact ansvver and the numerical approximation 
is 


0.000060 0.00012 
Ari -———- 


0.00054 0.00006 0.00964 1 

et nı na 7 vi, . (c) 
(A.3) 

Note that the final coefficient is only off by about one percent. VVith more 

numerical effort, this could be brought dov more. For the purposes of 

this article, it is important that the calculational pipeline vvorks, at least in 

principle. 


Appendix B. Alternative choice of numerators 


"The folloving choice of numerator factors for master integrals optimizes 
the povver series expansions of the coeflicients to some extent. "This choice 
is guided by the follovring principles: 


ə choose the numerators linear in the loop momenta to minimize the 
complexity of the Feynman parameter integrals: 


ə vvhen expressed in the master integral basis preferred by Reduze, the 
coefficient of the linear-numerator master integral comes vvith a nega- 
tive povver of € if possible, 


ə vvhen expressed in the master integral basis preferred by Reduze, the 
coefficient of the scalar version of the linear-numerator master integral 
comes vvith a, negative povver Of e, 


ə preserve permutation symmetry if possible. 
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This leads to the folloving list, vvhere the coordinates are in the basis of 


the integral topologies given in appendixlAppendix C 


INTİ130, (1,1,1,1,1,1,1,1,1,1,1,1,0,0, 0, 0, 0,0), pə - isl, 
INTEL 1.00, 0.00 Lab —İ: 
İNTES, İLL 1,1:1,11,1:1:T 10/00/0000 pə El, 
INTİ26, (1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0),pi - ql, 
INTİ25,11,1,1,1, 1:1,1:1/1,1, 1, 1/0,0,0)0,0,0k 0-ci “Bİ, 
INTİ23, (1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0), pə - isl, 
INTEL. ii in 0 ui Iİ, 
INTE7, 11/0/1111. 1,1 1, 10.0/0, 0/01 — İİ, 
INTİTL (1/0,1,1:1,1.1/10:1 1, 1:0/0, 10,0, -Eİ, 

İNTİB LUT 0.1... .0)0.0/0:0h:e—l) ıı 
INTA, 110,110, 11/1/1010 10,000, 0ksi bi, 

INTİ3,11, 1,0,1,1,1,1,0,1,1,0,1,0,1,1,0,1,0)p (4 — 2, 

INT 2. 10, 0:1,1,1,1,1.1:1/1:1/100.1 0/0 eb: 


ssl "This choice leads to a, significant reduction in the IBP computational vvork- 
Ço load, at the expense of introducing non-symmetric numerators. 


Appendix C. Basis 


— This appendix contains the main basis used in the reduction. The num- 

.- bering of the equations corresponds to the topologies in . ln each case, 

.——.—X the tvvelve entries given in the first line parametrize the tvvelve propagators 

P of the respective integral and the six entries in the second line the chosen 
numerators. Note that the numerators do not respect the symmetries of the 
diagram topology corresponding to the first tvvelve propagators, vvith the 
exception of topology 1. VVe have defined q — pı TL pə. 


He, 15, da, ls, 165, l4—ls, pi—ls, İ5 —l4, —ls-—q, —laq, —laq, —lsmq, 
5212757 o ə 700. 


la, İs 


q, İls-q, İn 


q, İs 


l4—l: 


Us, İs, da, İs, İ5—la, l6—İs, pi—le: İs 


İş — İs, İş — İs, la 


İs, IZ 


puls”, 0 


po), 


(C.2) 


İle: ls, dA, 5, 1a—la, da”, 16—İs: pi 


.. İq 


q, İls-rq, İs 


q, 4 


İs — pı, İs — pə, la 


Əl 


Do, İs 


Də, İs — İs, İs 


İşİ, 


15—l6 
(C.3) 


q, 


q, 
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He, ds, da, İs, pi-le, lala, 15—l6, 15—la, —ls-tlspə, —ls-q, —la-ıq, —lanq, 


pə, la — pə, ls — po), (C.4) 


İs — İş, İs — İs, la — İs, İs 


He, ds, da, ls, İs—la, 15-da, pi-ls, 15—l6, —ls—q, —lsmle-pə, —lsmq, —lsmla 


İs — pə, İs — İs, la — İs, la 


pı, la — pə, ls — po), (C.5) 


Hə, 15, la, İs, pi—ls, la—ls, la—la, p?—ls, —la 


le, İs — pə, la — po), (C.6) 


İs — İş, İs — İs, la — İs, İb 


İsq, 


q, —la-ıq, —la—lsq, —ls—ls-rq, 


q, —lals-ls, —laq, —ls—ls 


İs, 15, la, İs, pi—ls, l4—ls, la—ls, p?—ls, —lA 


İs — İs: İg — İş, is — İs, ə 


pə, la — pə, ls — po), (C-T) 


"q, 


İs, 15) la, İs, p?—ls, la—ls, la—lş, pi—ls, —la 


ls, la — pə, İs — po), (C.8) 


İs — İş, İs — İs, la — İs, İs 


"ə, 15, da, İs, la—ls, pi—ls, la—İls, —ls--q, —İş-bq, 


İs — 15, İs — pə, la — pə, İq 


(e.ə Pİ 16” pə), (C9) 


He, ds, da, ls, la—ls, pi-ls, l3—la, l3—lamls, 


l4 İ6 Fq, l5-l6—q, —— 


q, —lsq, —la—ls-q, la—ls-Tl6—pə, 


ls-rq, 


İsaq, —lsq, —l3—ls-tq, l3—ladlsls—q, 


pə, ls — pə, la — ls), (C.10) 


15— İs: le 10 İş, da 


l3—le-rq, İ3—la 


İ5-—l:—q, 


He, ds, da, ls, lalə, pi—ls, l4—İş, —İsq, —ls-rq, İsls—q, 


pə, 15 — le, ls — po), (C.11) 


la — İs, İs — pə, la — pr, 14 


He, 15, da, İs, p—le, pi-ls, l3—la, 1s—la, —İş 


İs — ls, ls — pə, İs — po, 14 


pı, la — pə, ls —pi), (C.12) 


q, —la lsrls, —ls—ls-q, —ls 


14—15—l6-q, 


İs-q, 


q, la ls-rq, İs İs-q, IZ 


He, ls, da, ls, la—la, p?—le, pi-la, l6—l5, —lə 


pı, 15 — pə, ls —pi), (C.13) 


İs — le, İş — pə, la — mə, 1: 


la—ls, —lsq, —la—lsmq, —l4—l6-iq, 


He, 15, da, İs, p?—le, pi—la, 13—la, 16—ls, —lə 


la — İş, la — ls, lA — ls, İs 


po, la — pə, ls — po), (C.14) 


q, —l4—ls-ls-Fpi, —l3—ls-bq, 


He, 15, la, İs, p?—ls, pi-la, 13—la, 15—l6, —lə 


İs — İs, İs — pə, la — pə: İs 


pı, ls — px, le —pi), (C.15) 
q, l3—larls—ls, —l3—lsmq, 


IZ isıq, 


4 İs q, 


He, ds, da, İs, p?—ls, pi-la, 16—ls, l3—la, —lə 


pə, 15— pn 16—pi), (C.16) 


İs — 16, İs — pə, la — İs, la 


ləl4—l:—ls”q, 


q, —ls—ls--q, —İstlA—ls-tq, 


Us, İs, la, İs, p—ls) 1s—la, pi—ls, 14—la, —la 


la — pə, la — ls: İ5 — pə, la 


pr, la — pə, 46 —pı), (C.17) 


lsl4—lsrls-pi, 


He, ds, la, İs, pi-ls, p?—ls, 13—la, l5—la, İs 


la—ls, —la-Fq, —ls”rl4”İaq, 
pə, ls — pə, le —pı), (C.18) 


İs — pə, İs — İs, la — İs, la 


öl 
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He, ls, da, İs, p?—ls, pi-ls, la—ls, lala, —ladlsmls, —lsq, —lslsmls, —l5—l6-q, 
İs — 16, 15 — le, la — pı, la — pə, la — pə, 15 — po), (C.19) 


İle, ls, la, İsi p2—le, la—ls, pi—ls) l5—la, —İs-tİs-rls, —la-q, —lsmrlals, —ls—le-kq, 
la — ls, ds — ls, la — pi, İs — pə, l4 — pə, İs — po), (C.20) 

İs, 15, la, İs, pi—la, la—la, pə—le, İserls, —la—lsmpi, —laq, —isLlempi, —l3—ls-pi, 
İs — pə, la — nə, İs — pə, 16 — pı, la — İs, İs —İs), (C.21) 


İle, ds, da, ls, İs—la, sls, p2—le, pi—la, —lsblspi, —ladq, —la—lstpi, —lslamls-ble, 
İs — İs, ls — pə, la — pə, İs — pi, ls — pə, l6 — pi), (C.22) 
Te, 15, da, İs, l6—pə, pi—la, 16—ls, İs—la, “İşıq, l3—la—ls, —ls-rls-rpi, —lsılsmpi, 
İş — pə, la — pə, İs — pı, ls — m, İs — pə, 14 —l6), (C.23) 
İls, 15, la, İs, p2—le, ls-le, l3—la, pi—la, —l3—lsmpi, —l4—ls-pi, —laq, —l3—ls—le-q, 
İs — pə, la — ls, la — pə, İs — pr, ls — pə) le — pi), (C.24) 
İs, ds, da, İs, pi—ls, İ5—la, pə2—ls, l3—la, —lamls-rls, —laq, —lsla—ls-pi, —lsıl4—ls—ls-rq, 
İs — İs, İs — İs, İs — ls, la — pi, la — pə, İ5 —o), (C.25) 


İlə, 15) la, İs, İs—la, l3—l4, pi—ls, po—ls, —lablsols, “laq, l3—la-bls—pə, 


— İş la — ls — İs “q, İs — İş, İs — ls, ls — İs, la — pı, la — pə, İs — po), 
(C.26) 


İle, 15, da, İs, İs-kls, po—ls, pi—la, İs—la, —la—ls”pi, —ladq, —lslamls-pə, “İslami, 
lA — ls, ls — pə, la — pə, İs — ni, ls — pə, l6 —pi), (C.27) 

İle, ls, la, İs, la—la, pi—la, p2—1s) İ5—ls, —ls-tq, 1s—la-rls—pə, —la—ls-kq, —İseklA—le-kpə, 
İs — pə, la — ls, l4 — ls, la — pə, İs — pz, le —pi), (C.28) 


Us, 15, la, İs, pi—la, la—la, pə—le, İs-bls, —la—lsmpi, —laq, —lsTl4—ls-po, —l3—ls—ls-rq, 
İs — pə, la — le, la — pz) ls — pi, ls — pə, le —pi), (C.29) 
İls, ls, la, İs, la—la, İsls, p2—ls, pi—la, —ls-q, —l4—ls-tni, —la-tlamls-pa, —ls-blA—ls-tpə, 
İs — pə, la — ls, la — pə, İs — pı, İs — pə, les — pi), (C.30) 
Us, 15, la, İs, İs—pr, polo, İs-tls, la—ls—pə, —l3—ls”pi, —ladq, larle—pə, —ls—ladq, 
İs — 16, la — pi, la — pə, İs — pi, ls — pə, $6 —pi), (C.31) 
İls, ls, la, ls, ls—pı, p2—ls, İsls, —ls—lspi, —labis-tpə, —latq, 14—15—ls, —ls—la-rq, 
İs — 16, la — pi, la — pə, İs — pı, İs — pə, le —pi), (C.32) 


He, 1s, la, İs, p?—le, 1s-ble, l3—pi, “laiq, ladls-pə, —l3—laıq, ladls-rls-pə, —l3—la—ls—lsq, 


la — ls, la — pı, la — pə, is — pi, ls — pə, le —pi), (C.33) 


ə2 


İle, İ5, la, İs, po—ls, l3—pi, İs-bls, la-kls—pə, —la—ls-pi, “laq, —l3—la-ıq, 


İs — le, la — pı, la — pə, is — pi, ls — po, 16 —pi). (C.34) 


Appendix D. List of ancillary files 


ə List of master integrals appearing in the form factor obtained by 
Reduze 


ə List of master integrals as suggested by Mint 


ə İBP reduction of all integrals appearing in integral topology 25 for the 
NN — 4 SYM form factor 
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